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TECHNICAL NOTE 3397 


AN EVALUATION OF NON-NEWTONIAN FLOW IN PIPE LINES 
By Rath, N. Weltmann 


SUMMARY 

An analysis is presented of a method for determining pressure losses 
due to the flow of non-Newtonian materials in pipe lines hy using "basic 
flow data obtained from measurements of flow curves, which are rate-of- 
shear - shearing-stress curves. The advantages of properly designed ro- 
tational viscometers over capillary viscometers for measuring these flow 
curves and interpreting them to obtain the basic flow parameters are dis- 
cussed. Dimensionless parameters are calculated from these basic flow- 
data and are used to construct a generalized friction diagram to describe 
the flow characteristics of Newtonian and non-Newtonian materials in pipe 
lines . 

Reported flow and pressure data are analyzed. Flow-pressure losses 
due to the flow of non-Newtonian materials in pipe lines, as calculated 
by the use of reported basic flow data and the generalized friction dia- 
gram, are compared with those obtained directly from reported pressure 
measurements in the pipe line. The in^jortance of considering the effect 
»; of tempera-ture on the flow parameters of non-Ne-wtonian materials, when 

making flow-pressure-loss calculations for -the flow in pipes, is men- 
tioned. The pipe transition pressure losses e:q)ected from non-Newtonian 
- materials are discussed. 


INTRODUCTION 

A great number of commercial materials, which heretofore have been 
considered to be Ne-wtonlan liquids, in reality show non-Ne-wtonian flow 
characteristics and thus have to be treated accordingly. Since many of 
these materials are passed through piping systems in their process of 
manufacture and application, the problem of the proper design of such a 
pipe-line system to carry non-Newtonian materials under given applied 
pressures and at desired flow rates has become of great Interest. 

Pipe-line resistance to the flow of Ne-wtonian liquids has been well 
understood for a long time and is frequently represented by a friction 
diagram, in -which a viscous friction factor is plotted as a function of 
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the Reynolds number for pipes of various surface roughness. Such diagrams 
can he found in the Literature (as, for instance, in ref. l) and extend 
over laminar and turbulent flow conditions. ^Hedstr5m (ref. 2) calculated 
a friction diagram for plastic materials of the Bingham type, based on a 
mathematical analysis of their flow curves (refs. 2 and 3) and expressed 
in temis of two basic flow parameters characterizing the nature of the 
flow. 


This paper presents a similar friction diagram, arranged for general 
use in calculating pressure losses in pipes induced by the resistance of 
Bingham plastics . This same friction diagr^ is extended to allow evalua- 
tion of pressure losses in pipe lines carrying pseudoplastic and dilatant 
materials; this diagram is based on an empirical analysis of the flow 
ciirves of -the materials and is expressed in te3mns of two empirical flow 
parameters characterizing the nature of the. flow. Experimental flow and 
pressure-loss data reported by various investigators for Bingham plastics 
(refs. 4 and 5) and for pseudoplastic materials (refs. 6 to 8) are used 
to check, the validity of this generalized friction diagram. 

An attempt is made to interpret the flow curves obtained for Bingham 
plastics and pseudoplastic and dilatant materials from measurements with 
a concentric-cylinder rotational viscometer and a capillary viscometer in 
order to obtain meaningful values for the flow parameters. A short dis- 
cussion of the flow characteristics of thixotropic materials is given to 
indicate the extent to which the generalized friction diagram may be used 
for those materials. The effect of temperature on the flow characteristics 
of all materials is mentioned. In addition,'” the pipe transition losses - 
such as occur in contractions, extensions, bends, orifices, and pipe en- 
trances, which are important for the proper design of a pipe-line system - 
are discussed for the case of non-Newtonian flow. 

The calculation of friction factors by the use of the characteristic 
flow parameters, which are obtained from noh -Newtonian viscometer flow 
curves, is demonstrated for a pipe line carrying a Bingham plastic and a 
pseudoplastic material at different flow rates. 


BASIC FLOW RELATIONS 
Newtonian Materials 

A Newtonian material is characterized by one material constant, 
namely, the Newtonian viscosity |j,. In laminar flow, the ratio of the 
shearing stress S to the rate of shear dv/dr is constant, and the 
viscosity jj. is defined by this constant relation so that 

- S 

^ dv/dr 


( 1 ) 
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(All symbols are defined in appendix A. ) Thus, the flow curve, which is 
the rate-of -shear - shearing-stress curve, is a straight line passing 
through the origin (fig. 1(a) ) . In turbulent flow, the shearing stress 
increases faster than the rate of shear (fig. 1(b)), and equation (l) is 
not valid. In fact, it is difficult to obtain the absolute viscosity of 
a Newtonian liquid from its turbulent -flow curve. 

The viscosity of a Newtonian material can be obtained from any "one- 
point" flow determination under laminar -flow conditions. 

If a concentric-cylinder rotational viscometer is used, equation (l) 
can be expressed so that 


S T/2rtR2h 



( 2 ) 


and thus. 


p = 


4 

4jthoi 


( 3 ) 


If a capillary viscometer or a pipe is used, equation (l) can be expressed 
(ref. 9) so that 


S _ ED/4L _ Sw 
dv/dr 8v/D 


and becomes Poiseuille's equation, where 

Pd2 

^ ” 32Lv 


( 4 ) 


( 5 ) 


The viscosity p has the units of poises if all the other variables are 
given in c.g.s. units. Equations (2) to (s) are derived under the assump- 
tion that slippage flow, end effects, and kinetic energy effects are 
negligible . 

Near the entrance of a capillary or pipe, the flow approaches the 
fully developed laminar -flow condition of equations (4) and (5) only 
asymptotically. In accordance with reference 10, this condition is 
achieved in laminar flow at an approximate entrance length of 
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Le = 0.029 Re[ii] D 


( 6 ) 


where Ee[|j.] = Dvp/n. 

For the description of pipe-line flow characteristics. It has heen 
customary to rewrite equation (5) in a dimensionless form, so as to re- 
late a friction factor (p to the Reynolds numher Re. For laminar flow 
of Newtonian fluids, then. 


where 

2PD 
“lp^ 

For turbulent flow of Newtonian liquids in pipes, empirical equations 
have heen established and are given (ref. l) as a function of the fric- 
tion factor (p. They are, for smooth pipes 


( 8 ) 


= 2 log (Re[pi3V^) - 0.8 


and for rough pipes 



+ 1.7 


( 9 ) 


( 10 ) 


Thus, the flow depends on the viscosity in smooth pipes but is in- 
dependent of viscosity in rough pipes under turbulent -flow conditions. 

In laminar flow the pressure drop is a function of the velocity (eq. (5)), 
while in turbulent flow it is a function of the square of the velocity 
(eqs. (9) and (lO) after substitution of the eaqiressions for Re[p] 
and <P ) . 


Bingham Plastics 


Bingham plastics are defined by Bingham* s equation (ref . 


ll), where 



(H) 


so that Bingham plastics have two characteristic material constants, 
namely, the plastic viscosity U and the yield value f. The flow curve 
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in laminar flow is a straight line with an intercept on the shearing- 
stress axis (fig. 2(a)). The yield-value intercept indicates that flow 
can start only if a shearing stress larger than the yield-value stress 
is applied. 

Experimental measurements at at least two points must be made under 
laminar-flow conditions to detennlne the plastic viscosity and yield 
value of a Bingham material. The instrument used for these determina- 
tions must be constructed in such a way as to make stress and shear cal- 
culations possible. 

In a concentric -cylinder rotational viscometer, equation (ll) can be 
integrated in accordance with Eeiner and Riwlln (ref. 12) so that 



When plotting flow curves, G is plotted against S. In the case of 
rotational viscometers, both these values are functions of R (eq. (2)). 
However, since the slope of the straight-line flow curves is independent 
of R, the viscosity values p. and U (eqs. (3) and (l2)) are not af- 
fected by the choice of R. This is different when c^culatlng the yield 
value from the stress intercept (Sq)|^ (fig- 2(b)), since (Sq^jj varies 

with the value of R (eq. (13)) that is used to calculate and 

for plotting the flow curves. The quantity T in equations (2) and (l2) 
is an applied torque and thus corresponds to a mean or Integrated value 
of Sj^. Mean values of rate of shear and shearing stress are approxi- 
mated by setting R = = (Rj^ + Rj,)/2. Thus, for convenience, these 

values will be called "mean rate of shear" and "mean shearing stress" and 
are used in plotting the flow curves. Then, and 

in equation (l2), and R = Eg^^ in equation (l5) . 

The shearing stress in a concentric-cylinder rotational viscometer 
(eq. ( 2 )) at a constant applied torque or mean shearing stress is larger 
at the Inside cylinder and decreases toward the outside cylinder. There- 
fore, the yield- value stress f is reached at the inside cy3J.nder first, 
and flow will start at this member when the local stress equals the yield 
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value Sf while the material will remaiu unsheared at the outside member. 
Hence, there is formed a "plug" whose radial thickness decreases with in- 
creasing applied stress and disappears when the local stress at the out- 
side cylinder equals the yield value. == 


Plug flow is always obtained with a concentric-cylinder viscometer 
at shearing stresses between R^f/R2 and R|f/R2 (fig. 2(b)), but the 

plug-flow portion of the flow curve can be held small by proper design 
of the viscometer. The angular velocity at which plug flow disappears 
is given as 


f 

^ =U 


_2Rg 



(14) 


Since equations (12) and (13) are valid only if the material is sheared 
over the complete radial thickness between R^j and R^, the plastic vis- 
cosity and yield value have to be determined from the straight portion of 
the flow curve (fig. 2(b)) at 

Only very few materleils are true Bingham bodies (fig. 2(b)). Most 
materials of this type will give a flow curve similar to the dashed curve 
in figure 2(b), which has a larger curved portion than would be produced 
by plug flow. The larger ciirved portion can be assumed in some cases to 
be caused by an elastic behavior of the materials, which is superimposed 
upon their viscous behavior. Such materials will nevertheless be called 
Bingham plastics as long as they can be approximated by a straight line 
so that their flow can be interpreted by the two Bingham flow parameters, 
namely, U and f . For these materials f might not be rigorously equal 
to, but rather a function of the yield value as defined by Bingham. Thus, 
to avoid misunderstanding and to stress its practical value rather than 
its physical meaning, the use of the term yield-value Intercept for f 
seems more appropriate for such materials. However, a meaningful choice 
of terms is possible only if the flow measurements permit a detailed 
analysis of the flow character of the material. 

Buckingham (ref. 13) integrated equation (U) for flow in a capil- 
lary viscometer and pipe, so that - 


U = 


RD^ 


32Lv 



(15) 


In deriving equations (12), (l3), and (15), slippage flow, end effects, 
and kinetic energy effects are again neglected. 
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The shearing stress in a pipe (eq. (4)), Sp = Pd/4L, varies along 

the pipe diameter d. At the wall in a pipe d = D and the shearing 
stress is = DP/4L, while in the center of the pipe Sp is always 
zero. The ratio 4If/DP = c is the ratio of the yield-value stress to 
the shearing stress at the pipe wall, and is equal to the ratio of plug 
diameter to pipe diameter. Since a plastic material always has a finite 
yield-value stress and since the shearing stress in the center of the 
pipe is always zero, plug flow disappears only at an infinite applied 
shearing stress and the flow curve hecomes a straight line only asymptot- 
ically (fig. 2(c)). For this reason, it is difficult to obtain U and 
f from capillary or pipe flow curvesj and to calculate them rigorously 
from equation (15) is even more difficult. Consequently, the two basic 
flow parameters U and f for Bingham type plastics are best obtained 
from rotational- viscometer flow curves. 

In analogy to Newtonian liquids, an equation for the friction factor 
for Bingham plastics flowing in pipe lines in laminar flow is 


<P = 



(16) 


if the plastic viscosity U is used, in the equation for the Reynolds 
number Re, so that Re[U] = Dvp/U. 


In accordance with an analysis by McMlllen (ref. 3) 
(ref. 20 j equation (16) can be rewritten so that 


and Eedstrom 


8P1 

^ “cRetUj 


64 

Re[U] 


PI 


PI 


" 8^ = ‘PoluJ ^ 


(17) 


where PI = fD/Uv is a dimensionless number, which will be called "plas- 
ticity number," since it Increases if the Bingham body is more highly 
plastic (larger yield value), decreases if the Bingham body becomes more 
viscous (lower yield value), and is zero for f = 0. Hedstrom (ref. 2) 
presents a graph which he determined by numerical calculation, giving 
<pRe[u] as function of PI. This graph has been used to extend the gen- 
eral friction diagram for the flow of Newtonian liquids in pipes to that 
of Bingham plastics. Such a generalized friction diagram for Newtonian 
liquids and Bingham plastics is presented in figure 3. In this diagram 
the friction factor <p is plotted as a function of the Reynolds ntimber 
Re, and the plasticity number is used as a parameter. 


Pseudoplastic and Dilatant Materials 

Flow curves of a pseudoplastic and a dllatant material are shown in 
figures 4(a) and (b), respectively. Power functions have been shown em- 
pirically to fit these flow curves and have been proposed as rheological 
equations for such materials (ref. 14), so that 
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dr 


(18) 


AlttLough this power function (eq. (l8)) Is not a true rheological 
equation, since it Is not derived from any physical concept, it can he 
regarded as an interpolation fonmila, and as such it is frequently very 
useful. In this equation, the quantity N indicates the structural 
change or the degree of pseudoplasticity and dllatancy and will he 
called "structure number." For Kewtonian liquids, W = 1; while W > 1 
for pseudoplastic, and If < 1 for dilatant materials. In many cases the 
value of If is constant over an appreciable range of rates of shear. 

The value of IT can then he determined from the slope of the straight- 
line plot of log (dv/dr) against log S. In cases where this plot is 

nonlinear, N = ^ and varies with the rate of shear and 

' d log S 

shearing stress. For ITewtonlan liquids where N = 1, equation (18) 
reverts to equation (l) and K 


1 ilquic 
= 1 /^. 


Substituting in equation (18) for S the expression of Sg for 
the rotational viscometer and that of Sp for "the capillary viscometer 

and pipe line, and assuming IT to he either constant over the range of 
rates of shearing obtained within the cross section of the viscometer 
test section for any one applied torque or pressure, or assuming N to rep- 
resent a mean value corresponding to the mean rate of shear, equation (18) 
becomes after integration (ref. 15} 


dr 



(19) 


where C has a different value for a rotational and a capillary viscom- 
eter. The rates of shear Gp and Gp (eqs. (2) and (4)) had been ob- 
tained under the assunption of constant viscosity and thus constant flow 
resistance. Since the flow resistance of non-Kewtonl^ materials varies 
as a function of shearing stress (eq. (18)) and the shearing stress 
varies within the confines of a rotational and capillary viscometer 
(eqs. ( 2 ) and (4)), the rates of shear produced within the sheared sam- 
ple of a non-Newtonian material will be different from those of a New- 
tonian liquid and will not only depend on the instrument used but also 
on the material being tested. The factor C in equation (l9), which 
depends on the material and the geometrical configuration of the viscom- 
eter, relates the rates of shear produced in non-Newtonian materials of 
the flow types as characterized by equation (10) to those calculated for 
Newtonian liquids at the corresponding shearing stress. Thus, the factor 
C will be called "rate -of- shear correlation factor." Then, for non- 
Newtonian materials having flow characteristics in accordance with equa- 
tion (is), the rates of shear are in a concentric-cylinder rota- 

tional viscometer for S = Sp, and Gjjj^p/Cp in a capillary viscometer 
or pipe for S = S-^. 
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The eagperimental ly ohtalned K, if it is not constant hut is a 
function of rate of shear, represents a mean value corresponding to the 
mean rate of shear only, if its change over the sheared width of the 
sample is very small. In a rotational viscometer of proper design, the 
change in shearing stress and thus of rate of shear within the annula 
between hoh and cup is small, so that, even for large values of dU/^, 
the change in H over (R^ - will he small. Thus N, wh^. experi- 
mentally obtained with a rotational viscometer, mif^t he considered a 
mean value at the measured flow condition, even if it changes substan- 
tially with a change in flow conditions, and equation (l9) is valid. 
However, in a capillary viscometer or pipe line the change in shearing 
stress along the diameter is always large, changing from zero to S^, so 

that the change of H over D is small only if dS/dS is small. Thus 
H, when experimentally obtained with a capillary viscometer, is often not 
a mean value, and equation (19) is Hien only an approximation. In this 
case the experiment ally obtained values for II and Cp are not any more 

constants of the material, alone but will depend on the dimensions of the 
viscometer. 

The quantity K in equations (iS) and (19) is a constant for a 
given material, but it has no physical meaning and its dimensions change 
as a function of N and thus afe different for each material. To give 
K physical meaning, equations (l8) and (19) might be rewritten by mahing 
use of such established flow constants as the relaxation time t and the 
shear modulus M, so that 


^ = G ^ 1 /S f ( 20 ) 

dr C T \M/ ^ 

It would be best to use basic material constants such as t and M for 
the flow parameters in the analytic process of establishing an equation 
for the friction factor as a function of Reynolds number for pseudoplastic 
and dilatant materials. However, since equation (20) has no theoretical 
or experimental basis, and since H is not expressed as a function of 
established material constants, this equation cannot be used to inter- 
pret the flow curve in terms of t and M, even if these two material 
constants are obtained from independent measurements. Thus, for the 
present, it seems useful and ei^edient to employ an empirical flow pa- 
rameter, namely, the apparent viscosity tj. The apparent viscosity q 
is defined as the ratio of shearing stress, to rate of shear, so that 

n = f- (21) 

For a concentric-cylinder rotational viscometer. 
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n = 



4irhco 


C 


R 


where 


f 1 - (Rb/Rc)2^ 1 

'/p\2(N-l) 

|N[1 - (Rb/Rj2jj 

W 


( 22 ) 


(23) 


and approaches 1 for Rb/^c ** capillary visconffiter or pipe, 

1 = ^ ( 24 ) 

where 

<=P = WT3 (2=) 

Thus, two empirical flow parameters can he stipulated for pseudoplastic 
and dilatant materials, namely the apparent viscosity t) (eq. (2l) ) and 
the structure nuniber N (eq. (l9)). However, for non-Newtonian materials, 
q is not a constant, but varies continuously with the applied rate of 
shear (eq. (21) ) . Frequently N also varies with the rate of shear, but 
from experiments it was found that in many of these cases N can be ap- 
proximated by a different -constant value over different rate-of -shear 
interveils. Nevertheless, both empirical flow parameters q and N are 
meaningful only if they are obtained at the same local rate of shear at 
which the flow resistance is to be determined. These "local flow param- 
eters" must be obtained from a flow curve extending over a wide range of 
rates of shear. An .approximate local apparent viscosity might be ob- 
tained for any specified rate of shear by extrapolation, if it is outside 
the flow- curve measurement, by using a logarithmic plot of equation (19 ) 
to determine N at the highest measured rates of shear. Assuming then 
N to be constant for still higher rates of shear, S can be extrapolated 
for any specified rate of shear by using the same plot. Then the local 
apparent viscosity is calcualted from equation (2l). 

In a capillary viscometer and in a pipe line, the mean rate of shear 
^ j/Cp is equal to _ 2v(N + 3)/D, (eqs. (4) and (25)). The mean rate of 

shear in a concentric-cylinder rotational viscometer is 

(eqs. (2) and (23) for R = R^^). 




u 
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A rotational viscometer is the preferred instorument for measuring 
the flow curves of pseudoplastic and dilatant materials, since, as pre- 
viously analyzed, the experimentally obtained If represents in most 
cases a mean value, and also because the rotational-viscometer mean-rate- 
of- shear coirelatlon factor Cr (at E = R^y) vlll in most practical 

cases approach 1 very closely (see fig. 5), so that the apparent viscos- 
ity (eq. ( 22 )) and the flow curve (eq. (19)) are almost independent of 
H. In a pipe line or caplU.ary viscometer, however, the experimentally 
obtained R often does not represent a mean value, and the rate-of-shear 
correlation factor Cp can differ from 1 considerably, as is shown in 

figure 5. In fact, the factor Cp is very dependent on N, so that the 

apparent viscosity (eq. (24)) and the flow curve (eq. (19)) depend greatly 
upon the validity of the determination of F. The apparent viscosity as 
given in equation (24) is the apparent viscosity only if the If substi- 
tuted in equation (24) is a mean value. Otherwise, effective apparent 
viscosities are obtained, which depend on the dimensions of the viscom- 
eter. Another disadvantage of using the capillary viscometer for the 
determination of the flow character of non-Newtonian materials is that 
the flow curves of a pseudoplastic material (fig. 4(a)) and of a Bin^am 
plastic (fig. 2(c)), irfien plotting Gr^^p against S^, are so much alike 

as to make it difficult to differentiate between one and the other flow 
type from a study of such flow curves, and thus to obtain meaningful mean 
rates of shear and flow parameters for use in scaling the measured re- 
sults to the flow conditions required for any specific application. 

In analogy to Newtonian liquids, an equation for the friction factor 
<p for pseudoplastic and dilatant materials for laminar flow in pipes is 

64 /n + 3N r /N + 3\ 

if the local apparent • viscosity q = SCp/G^^p (eq. (2l)) is used in the 
equation for the Reynolds number so that Re[q] = Dvp/q. 

In the case that the apparent viscosity is measured in a capillary 
viscometer, q = measured value of N might not rep- 

resent a mean value, so that Cq can differ from Cp. Do obtain the 
local apparent viscosity, an approximate mean rate of shear Gjjj^p/Cq 
I s substituted for Cm,c/Cc^ "the approximate local apparent viscosity 
is q = S^Cq/Gjjj^^P where Gj^^p= 8v/d. Then the equation for the fric- 
tion factor (eq. (26)) can be rewritten to read <p =* <p[qs]G^c/^P 
qg is defined as qg = S/Gjjj^p. Thus, for Cq = Cp, equation (26) 
reduces to 
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In this equation, Tig is not a viscosity, and therefore Re[T]g] is not 

a Reynolds number in accordance with its definition, hfut equation (27) 
can he useful. However, the condition that Cq = Cp is fulfilled only 

if the measured H represents a mean value, which it does only under 
the conditions that dN/dS is small. In other cases, Cq might 

differ considerably from Cp, and thus the use of a capillary viscometer 
for the measurement of ti and <p is questionable. 

In the case that t) is measured in a rotational viscometer, 

^ ~ ®m,R*-R/^m,R-> value of H, when measured in a properly de- 

signed instrument, most often represents a mean value. Therefore, the 
mean value of Cp for the pipe can he calculated (fig. 5). To obtain 

the local apparent viscosity, the mean pipe rate of shear <^^p/Cp is 
substituted for G^^p/Cp, and the local apparent viscosity is 
^ “ ®m,R*^p/*^m,P* equation for the friction factor reduces then to 

equation (27). 

Equation (26) is used to extend the friction diagram (fig. 3) to 
the flow of pseudoplastic and dilatant materials by using the structure 
number H as parameter, as indicated by the dashed lines. Since equa- 
tion (26) was derived from an enipirlcal formula and not from a rheo- 
logical equation, such a diagram should be considered only as a useful 
approximate method for calculating flow pressure losses of such materials 
in pipe lines. 


Thixotropic Materials 

Flow curves of a plastic, a pseudoplastic, and a dilatant material 
of thixotropic character are shown in figures 6(a), (b), and (c), respec- 
tively. Thixotropic materials change in flow resistance with rate of 
shear and shear duration under isothermal conditions. Flow curves of 
thixotropic materials are characterized by a hysteresis loop, when the 
flow curve is obtained with a viscometer that is capable of applying to 
the same sample an increasing shear rate followed by a decreasing shear 
rate. In the case of thixotropic behavior, the "down curve" obtained at 
decreasing rates of shear win be shifted towards smaller stresses com- 
pared with those of the "up curve" determined at increasing rates of 
shear. This Indicates that the flow resistance decreases with increasing 
rates of- shear and shear duration. Upon rest, the flow curve will shift 
back, which indicates that the flow resistance increases to its original 
value when the material is left unsheetred. 
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During the measurement of the up curve, the structure of the mate- 
rial is continuously "broken down, thus changing the flow character of the 
material, so that each point on the up curve mi^t "be considered to rep- 
resent a point on a different flow curve. According to the literature 
(ref. 16), each down curve, when measured rapidly, is said to represent 
a flow curve of the material. Since the down cxrrve shares one point, 
namely the top point, with the ttp curve, the down-flow curve describes 
the flow behavior that characterizes the material at this point of its 
shear history. It has been found (ref. 16) that the position of this 
point and thus that of the down curve depends on two factors, namely, 
the value of the highest rekte of shear applied to the sample, which 
will be called "top rate of shear" of the flow curve, and the time of 
shear duration before the down cttrve is commenced. 

Most capillary viscometers are incapable of applying more than one 
rate of shear to the same sample^ therefore, they do not produce hyster- 
esis loops and are unsulted for the measurement of thixotropic materials. 
Rotational viscometers, however, are capable of detecting hysteresis in 
flow curves and thus are well stilted for the measurement of thixotropic 
materials . 

The flow parameters are obtained from the down-flow curves in a 
similar manner as was done for materials of non-thixotropic character. 

In the case of thixotropic materials where the down curve characterizes 
the material as a Bln^am plastic (fig. 6(a)), the flow parameters are 
the plastic viscosity U and the j^eld-value intercept f, and equa- 
tions ( 12 ), (13)^ ( 17 ) applicable. Da the case of thixotropic 

materials where the down curve characterizes the materials as being 
pseudoplastic or d1 latent (figs. 6(b) and (c), respectively), the flow 
parameters are the local apparent viscosity t) and the local structure 
number N, and equations (I 8 ), (2l), and (26) are applicable. However, 
for a1 1 thixotropic materials, the flow parameters obtained from a down- 
flow curve describe the flow condition of the material only if it has 
been subjected to a specific shear treatment that is identical to that 
which the sample experienced in the viscometer before the down-curve 
measurement. For this reason, the flow parameters thus obtained will be 
called "specific flow parameters." For a meaningful interpretation they 
should always be marked by the top rate of shear and the shear duration 
to which the material weis subjected before the measurement. 

For thixotropic plastic materials, an empirical extrapolation method 
has been suggested (ref. 16), which permits the calculation of the spe- 
cific flow parameters U and f for any specified value of top rate of 
shear and shear duration. To do this, three flow curves have to be ob- 
tained, two at two values of top rate of shear, and a third at one of 
those top rates of shear but after a d if ferent time of shear duration. 
Such an extrapolation method permits the calculation of the specific flow 
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parameters for shear conditions, vhich might not easily he obtained with 
viscometers. For pseudoplastic and dllatant materials of thixotropic 
character, no extrapolation method has been established. Thus, it is 
necessary to make the measurement -of the flow curve under the Identical 
specified conditions of top rate of shear and shear duration to which 
the material is subjected in the pipe line. Only then is it possible to 
obtain meaningful flow-pressure-loss calculations from the interpretation 
of such a flow curve. 

The use of such specific flow parameters to obtain flow pressure 
losses in pipe lines by employing the friction diagram can best be con- 
sidered a rnugh approximation method. The flow of thixotropic materials 
in pipe lines is very con^ilex. In any pipe line the shearing stress 
changes from zero at the center of the pipe line to a m aximum at the pipe 
wall, with coirespondlngly large variations in rate of shear and shear 
duration along the radius. Since the specific flow parameters of most 
thixotropic materials are greatly affected by“ changes in rate of shear 
and shear duration, they will vary substantially along the entire cross 
section. Thus, in the case of thixotropic materials, it is rather dif- 
ficult from pipe-line flow-rate measurements to plot a curve of mean 
rate of shear as a function of mean shearing stress, and to determine 
meaningful values of mean rate of shear and of mean shear duration at 
which to measure the specific flow parameters.. For this reason it seems 
that experimental data would be useful to establish mean values of shear 
conditions as a function of flow rate at which to calculate the specific 
flow parameters for determining flow-pressxire "losses of thixotropic ma- 
terials in pipe lines by the use of the friction diagram. 


GEHERALIZKD FRICTION DIAGRAM 
Newtonian Materials 

For Newtonian liquids the plasticity number is zero. The lowest 
solid curve in the friction diagram (fig. 3) is for smooth pipes. At 
Re[pl <2000, the flow is laminar and equation (?) is valid. For 
turbulent flow, equation (9) is valid for smooth pipes and equation (lO) 
for pipes with an extremely rough surface of D/k < 20. For flow in 
pipes of Intermediate surface roughness, lines between those two curves 
can be plotted by using equation (lO) or can be taken from the litera- 
ture (ref. l). 


Bingham Plastics 

For Bingham plastics, the lines of different plasticity numbers 
(fig. 3) are representative of the laminar-flow region. As long as the 
plasticity number and Reynolds number can be calculated, the friction 
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factor can be obtained from the diagram (fig. 3). An evaluation of ex- 
perimental data reported for smooth pipe lines (refs. 4 and 5) indicates 
that the flow of a Bingham plastic in a pipe line becomes turbulent at a 
Reynolds number at which the line of respective plasticity number inter- 
sects the turbulent branches of the curves for Ifewtonian flow. For in- 
stance, at a plasticity number PI of 100, turbulence will start at 
Re[u] > 4.7x10- in smooth pipes and accordingly would be expected to 
stant at Re[Uj > 1.4x10^ in the extremely rough pipes. The transition 
Reynolds number for Bingham plastics increases with increasing plasticity 
numbers . 


Pseudoplastic and Dilatant Materials 

For pseudoplastic and dilatant materials, the lines of different 
structure number (fig. 3) are representative of laminar flow. An approx- 
imate friction factor can be obtained if the local structure number and 
the local apparent viscosity (eq. (2l)) can be calculated. In analogy 
to plastic flow, it is assumed that turbulent flow for such materials 
starts at the Reynolds number Re[Tj] at which the line of respective 
st 2 mcture number intercepts or meets the turbulent branches of the curves 
for Ifewtonian flow. Experiment ai data (refs. 7 and 8) indicate the va- 
lidity of such an assunption. The friction diagram for pseudoplastic and 
dilatant materials is similar to that for Bin^am plastics, except that 
different flow parameters are used in its interpretation and that the 
spacing of the lines is different. Since the spacing of the lines is 
uniform throughout, all friction factors can be obtained Just as readily 
(eq. (26)) by using the lines for the Newtonian friction factor to ob- 
tain <Pq[ti] and figure 5 to obtain Cp. However, the onset of turbu- 
lence can be determined only by using the lines of respective structure 
number in the friction diagram. 


CALCULATIOB OF PRESSURE LOSSES IN PIPE LINES 
Uniform Pipe Line 

For designing pipe-line systems, the flow rate and applied pressure 
are frequently specified, so that the pipe lines have to be dimensioned 
to produce a min:fTnuTn pressizre loss for otherwise optimum operational con- 
ditions. The pressure loss in "uniform pipe lines" - that is, in pipe 
lines of constant cross-sectional area - is due entirely to the flow 
resistance of the materials carried in that pipe. This "flow-pressure 
loss" can be readily calculated -when the friction factor is known (eqs. 
(7), (17), and (26)). 

The relation between the friction factor and the flow-pressure loss 
is the same for materials of any flow character by the definition of the 
friction factor (eq. (8)). The equation for the flow-pressure loss (ref. 
1) is 
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APp. = 



(L - Lg) 
D 


tp 


(28) 


Equation (6) for the entrance length Lg is strictly valid only for 

laminar flow of Nefwtonian materials. Under all other flow conditions^ 
equation (6) is at hest an approximation, hut useful enough for most prac- 
tical pipe-line considerations. 


Konuniform Pipe Line 

Practical pipe-line systems consist hy necessity of uniform and non- 
uniform pipe-line portions. Nonuniformity can he produced in the pipe 
line hy various means, such as contractions, extensions, hends. Inserts, 
orifices, entrances, and others. The pressure loss in the nonuniform 
portions of the pipe line, which will he called "pipe transitions," is 
again due to the flow resistance of the material passing through the non- 
uniformity. However, this "pipe transition loss" cannot he readily cal- 
culated, since the flow is nonuniform, so that a very complex flow pat- 
teni results. The pipe transition loss is 

APj, = p ^ Cl (29) 

For Newtonian liquids the pressure-loss coefficients Cl have heen 

determined e 2 q)erlmentally for various types of pipe transitions and can 
he found in the literature (refs. 1 and 10). For convenience, some 
values of Cl for Newtonian liquids are given in table I. 

The pressvire losses are additive, so that the total pressure loss in 
a pipe -line system is 


v^ r (^ ” Le) /• . 

APj. + = p I^Cl + • p <Pj (30) 

In practical pipe-line systems of large L/D, when (L - /D > > 
the pipe transition loss can frequently he neglected. 

From an analysis of the generalized friction diagram (fig. 3), it 
seems that the pipe transition losses for the flow of non-Newtonian ma- 
terials could become rather large, since the friction factor in laminar 
flow can he a multiple of that for a Newtonian liquid of the same vis- 
cosity. In turbulent flow, however, the friction factors would he equal 
to those of Newtonian liquids. Thus, it would seem that the pipe transi- 
tion losses of non-Newhonian materials will greatly depend upon whether 
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laminar or tiorbulent flow prevails In the pipe transitions. In addition, 
non-Newtonian materials such as Bingham plastics require that the shear- 
ing stress exceed the yield-value stress before flow can- start. This 
requirement tends to make the pipe transition losses not only dependent 
upon the configuration of the pipe transition but also upon the magnitude 
of the yield-value stress, since plug flow will occur whenever the shear- 
ing stress decreapes below the yield-value stress, and an Increase in the 
amount of plug flow will Increase the pressure losses as long as laminar 
flow prevails. Since dilatant materials (fig. 4(b)) Increase in flow 
resistance with increasing rates of shear, any pipe transition that con- 
tains a constriction will produce a large pipe transition loss, because 
the material is subjected to a relatively high shear in such a 
constriction. 

Only one reference was found (ref. 3) presenting pipe transition 
losses of non-Newtonian materials. In this particular reference, the 
flow of a Bln^am plastic was studied. The pipe transition losses re- 
ported were for a pipe contraction and were many times larger than would 
have been expected for a Newtonian material of the same viscosity sub- 
jected to the same flow conditions. This trend seems in accord with the 
preceding analysis. 


EVAIUATION OF EXEEBIMENTAL DATA IN EELATION TO FIGI3BE 3 

Bingham ELastics 

Experimental data of the flow-pressure losses of Bingham plastics 
in smooth pipe lines are reported in references 4 and 5. They are used 
to check figure 3. HedstrSm (ref. 2) has used some of the pressure-drop 
data reported in reference 4 to establish the onset of turbulence for 
Bingham plastics flowing in pipes. For this reason HedstrSm has computed 
the plastic viscosity and yield value of the material used in this refer- 
ence, which is a cement rock suspension. The same values of U and f 
were used to obtain the data in figure 7. The other investigator (ref. 

5^ studied clay in water suspensions at different solid contents. He 
measured the pressure loss in the pipe line over a length of 20 feet and 
obtained the two basic flow parameters, U and f, at the same tempera- 
ture from flow curves produced with a rotational viscometer. The data 
used in figure 7 Include a range of 0.08 S U S 0.235 poise, 0 < f S 276 
dynes per square centimeter, 1.91 < D <7.60 centimeters, 0.01 < v <690 
centimeters per second, and 1.01 < p <1.52 grams per cubic centimeter, 
giving 0 < PI £.480,000 and 0.74 < Re[U] < 161,000. The friction factor 
^experimental from the experimental pressure-drop data measured 

in the pipe lines is con^ared in figure 7 with the friction factor 
'^diagram'' calculated by use of the generalized friction diagram 

and the two basic flow parameters. The agreement of both friction factors 
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up to cp » 10^ shows that the generalized friction diagram is valid for 
evaluating the flow in pipe lines of Bingham plastics from their two 
basic flow parameters. 

Above tp » la^ the plasticity number was high (Pl > 600) and the 
flow velocity was in most Instances low. In fact; most of the points 
which are more than 50 percent off were taken at a flow velocity of 
V <0.3 centimeter per second. Since, as was mentioned before, only 
very few materials are true Bingham bodies but most show an elastic flow 
behavior (dashed curve in fig. 2(b)) at the low shearing stresses, their 
yield-value intercept and plastic viscosity are constant only above a 
certain shearing stress S** which is Isrger than the calculated yield- 
value stress. Since it will be shown later that for values of PI > 600 
the accTiracy of the friction factor is independent of that of the plastic 
viscosity, only the change in yield-value intercept at the low velocities 
is considered. 


At shearing stresses between S* and S** the yield-value inter- 
cept increases continuously with increasing rate of shear until it reaches 
at S** and (dv/dr)^ its maximum constant value of f. Thus, at flow 
velocities corresponding to rates of shear of less than (dv/dr)** the 
flow parameter f is larger than the local yield-value intercept at the 
given local rate of shear. Consequently, the calculated plasticity num- 
ber, being proportional to f, is also too large and so is the friction 
factor obtained from figure 3. More accurate ffiction factors 

can be obtained from the friction diagram even for low flow velocities 
corresponding to a rate of shear of less than (dv/dr)** if the plasticity 
number is calculated by using the local yield-value intercept, namely, 
that which is characteristic for the rate of shear prevailing in the pipe 
line. The-local yield-value Intercept, which should replace f in the 
calculation in such cases, can be approximated by using the Intercept on 
the stress axis, which is obtained by constructing a tangent on the flow 
curve at the point of local rate of shear. 


An analysis of the accuracy of the friction-factor determination as 
a function of the accuracy of the measurements of U and f, and of flow 
velocity, density, and pipe dimensions, can be made with the assistance 
of figure 8. Figure 8 is a logarithmic plot pf the plasticity number PI 
against «p/<PofUj» At plasticity numbers of PI > 600, this logarithmic 

Aq>o [U j 


^ Atp 1 /api\ . 

Am API A9 o[u3 

= — — + After substitution of the ex- 

(p PI TqM 

pressions for PI and <Pq[u], it is found that the accuracy of the fric- 


relation is linear, so 


a slope of a = 1.0, 


tion factor depends on the accuracy of meas-urement of the yield value, 
flow velocity, and density and that it is independent of the accuracy of 
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measurement of the plastic viscosity and pipe dimensions. Then^ the 
accuracy of the flow-pressure-loss calculations (eq. (28)) at PI ^ 600 
depends entirely on the accuracy of the yield-value determinationj if 
the errors in the pipe-line dimensions and the pipe entrance length are 
assumed negligible. 

In the range of 20 < PI < 600, another linear approximation can be 
made (fig. 8). Since this line has a slope of ex « 1.2, it is found 
after substitution that the accuracy of the friction factor in this re- 
gion of plasticity numbers depends on the accuracy of measurement of both 
flow parameters and of density, flow velocity, and pipe diameter. How- 
ever, considering the Influence of the flow parameters alone, it is no- 
ticed that the accuracy of the yield-value determination has to be about 
five times better than that of the plastic viscosity measurement to pro- 
duce an equal error in the determination of the friction factor. With 
decreasing plasticity numbers PI < 20, the accuracy of the friction 
factor depends less and less on that of the yield value but to an in- 
creasing extent on that of the plastic viscosity measurement. Both flow 
parameters have to be determined to an equal accuracy for an equal error 
in the friction factor at a plasticity number of about 5, since at this 
point the tangent on the curve in figure 8 has a slope of tx » 2.0. 

With f\irther decreasing plasticity numbers, the influence of the yield- 
val-ue determination on the accuracy of the friction factor decreases 
rapidly and is insignificant at PI 3? 1.0. 


Pseudoplastic Materials 

A few papers in the literature (refs. 6 to 8) deal with suggestions 
for obtaining the pressure loss due to the flow of pseudoplastic materials 
in pipe lines from flow data obtained in a viscameter. In laminar and tur- 
bulent flow, good correlations have been obtained between the pressinre loss 
directly measured and the one calculated by using the Newtonian friction 
diagram after substituting a local apparent viscosity into the ej^resslon 
for the Re 3 raolds number, where the local apparent viscosity was calculated 
in accordance with equation (18). 

In laminar flow, the use of only the one flow parameter of local 
apparent viscosity assumes that the term (N + 3)/4 in equation (26) is 
close to 1 or that this term was cancel-led out by the use of an appro- 
priate equation for so that equation (27) was applicable. The former 
is valid only for structure numbers of N approaching 1. The latter is 
valid whenever rotational-viscometer flow curves are used to calculate 
the local apparent viscosity at a given pipe rate of shear and is approx- 
imated if the calculation is made from capillary- vis cometer measurements. 
Since the reported pressure-loss correlations were good, even in laminar 
flow, it is assumed that either one or both of the two conditions were 
satisfied in the eaqjeriments . 



20 


NACA TN 3397 


For Infinitely high, rates of shearj approximate values of local 
apparent viscosity can also he obtained hy exhrapolation from a plot of 
inverse rate of shear against local apparent viscosity. By assuming 
this plot to he linear, extrapolation to infinite rate of shear yields 
an intercept on the apparent- viscosity axis equal to the local apparent 
viscosity at infinite rate of shear. In turbulent flow, when the rates 
of shear are so high as to approximate the flow conditions at infinite 
rate of shear, such a ’'residual apparent viscosity" can he used instead 
of the local apparent viscosity to calculate Re[Tij. This has the ad- 
vantage that the knowledge of the local mean rate of shear in the pipe 
is not required. The residual apparent viscosity has been used in ref- 
erences 7 and 8 in the turbulent -flow evaluations. 

The accuracy of the friction factor for pseudoplastic and dilatant 
materials as a function of the accuracy of the measurements of the flow 
parameters q and H, and of the flow velocity, density, and pipe dimen- 
sions can be determined by using equation (26). After substitution of 
the expression for <PQ[n] equation (26), it is found that the accu- 
racy of the friction factor depends on the accuracy of measurements of 
both flow parameters and of density, flow velocity, and pipe diameter. 
However, considering the influence of the two flow parameters alone, it 
is noticed that the acc\uracy of the detemiination of the local struct\u*e 
number JJ can be less than that of the local apparent viscosity q by 
a factor of l/[l + (3/K)] to produce an equal error in the deteirmina- 
tion of the friction factor. Thus, for N » 2, the measurement of the 
structure number can be about 2.5 times less accurate than the apparent- 
viscosity determination to result in the same error in the friction fac- 
tor and flow-pressTire-loss calculations. 


temperature combiderations 

The flow parameters are functions of temperature. Therefore, the 
flow parauieters used for flow-pressure-loss calculations should prefer- 
ably be measured under the same temperature conditions as prevail in the 
pipe line where the pressure loss occurs. Frequently, this temperature 
cannot be easily determined, since it might depend on environmental fac- 
tors and thus comprise a complete range of temperatures, rather than a 
single value. In that case it is necessary 'to measure the flow param- 
eters as a function of temperature. Frequently, measurements of the flow 
parameters at two temperatures will permit interpolation or extrapolation 
over the. required temperature range with sufficient accuracy to allow 
calculation of the corresponding friction factors and flow-pressure 
losses. In Newtonian flow of liquids, temperature Increases will, in 
most cases, cause a decrease in viscosity and thus will effect a decrease 
in pressure loss in the pipe line. However, in non-Newtonian flow, an 
increase in pressure loss can be expected more frequently, since temper- 
ature Increases can affect either one or both of the flow parameters in 
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either direction. They can cause, for instance, an Increase in yield 
value in the case of a Bingham plastic, and thus effect an increase in 
pressure loss, in spite of a decrease in viscosity. 


CONCECJDING EEMARKS 

The importance of estahlishing Basic flow parameters, which for non- 
Newtonian materials can he obtained only from flow curves (rate-of -shear - 
shearing-stress curves), has Been discussed. Use has Been made of these 
characteristic flow data for calculating dimensionless parameters, which 
are presented in a generalized friction diagram. The purpose of this dia- 
gram is to assist in the evaluation of pressure losses resulting from 
the flow of non-Newtonian materials in pipe lines. 

A concentric -cylinder rotational viscometer of proper design has 
Been shown to have advantages over the capillary viscometer when measur- 
ing the flow of non-Newtonian materials to obtain the flow parameters. 

In fact, in some cases, the capillaiy- vis cometer flow curves do not per- 
mit a valid interpretation of the flow character of the materials. 

The two Basic flow parameters for Bingham plastics are the plastic 
viscosity and the yield value. The local apparent viscosity and the 
structure nttmber are used as characteristic flow parameters for pseudo- 
plastic and dllatant materials. For thixotropic materials, a method of 
using specific flow parameters corresponding to the type of flow curve 
and to the nature of the material is suggested. 

In the case of Bingham plastics, it has Been shown that for highly 
plastic materials the accuracy of the yield-value measurements alone 
determines the accuracy of a pressure-loss calculation when made By use 
of the friction diagram, while the accuracy of measurement of the plastic 
viscosity and the flow velocity Become in^jortant only for less plastic 
materials. In the case of pseudoplastic and dllatant materials, it has 
Been shown that the accuracy of a pressure-loss calculation, when made By 
use of the friction diagram, depends on the accuracy of measurement of 
the local apparent viscosity and the flow velocity and only to a lesser 
degree on the accuracy of the determination of the local structiire numBer. 

The validity of the generalized friction diagram for descriBing flow 
of Bingham plastics in pipe lines is shown By an evaluation of experi- 
mental data of two Independent investigators. Reference is made to other 
reported experimental data, which indicate that the generalized friction 
diagram is valid also for evaluating the flow of pseudoplastic materials 
in pipe lines. 

Considerations of the effect of temperature on the various flow pa- 
rameters suggest that temperature can affect pressure losses in pipe 
lines carrying non-Newtonian materials more adversely than in the case of 
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Wewtori±an flow. IThus, non-Newtonian flow pafameterB must te determined 
at various temperatures to establish their trend as a function of temper- 
atiore . The flow parameters corresponding to_ the approximate range of 
temperatures prevailing in the pipe line where the pressure loss occurs" 
can then he selected for the calculations of "the friction factors and 
flow-pressure losses. 

An analysis of the basic flow data in reference to the friction fac- 
tors suggests that pipe transition losses of non -Newtonian materials will 
in most cases, especially where the flow is laminar in the pipe transi- 
tion, he larger than those expected for Newtonian liquids of the same 
viscosity and under identical flow conditions. For this reason, transi- 
tion pressure-loss coefficients for the flow of non-Newtonian materials 
in pipe transitions are required. Since th^ cannot he obtained by cal- 
culation, they have to he established by experiment. 

Numerical examples of friction-factor calculations by the use of 
flow curves and the generalized friction diagram are given in appendix B 
for a Bingham plastic and a pseudoplastic material flowing in a pipe line 
at different flow, rates. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 26, 1954. 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report. The dimensions of 
all dimensional quantities are indicated. 

C rate-of -shear correlation factor 

Cq contraction coefficient 

pressure-loss coefficient in a pipe transition 

c ratio of plug to pipe diameter 

D diameter of pipe or capillary, [L] 

Dq orifice diameter, orifice slot width, [Lj 

d any diameter of pipe or capillary, 0<d<D, [l] 

dv/dr rate of shear or velocity gradient, [T“^J 
f yield value or yield- value intercept, [mL"1T"2J 

G rate of shear, assuming constant viscosity, [T“^] 

h hei^t of inner cylinder (hoh), [l] 

K proportionality constant 

k grain diameter of surface roughness in pipe, [L ] 

L length of capillary or pipe, [Lj 

M shear modulus, 

N structure nimiber 

P, AP pressure difference over given length L, 

AP_ flow-pressure loss, 

r 

AP^ pressure loss in pipe transition, 
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PI plasticity number, PI = fD/Uv 

R any radius, of rotational viscometer, Rb<R<Rc^ [L] 

Rav average radius between R-^^ and R^, R^^^ = -£-- ^ [L ] 

Rt, radius of inner cylinder (bob), [L ] 

Rq radius of outer cylinder (cup), [L] 

Re Reynolds number Re = Dvp/viscosity 

r radius of bend of pipe, [l] 

S shearing stress, 

Sq stress intercept at G = 0, 

T torque, 

Tq torque intercept at cu = 0, 

U plastic viscosity, 

« 

V mean velocity of flow (calculated from flow rate), [lT"^] 

a slope of lines in fig. 8 

Ti apparent viscosity, tj = s/(dv/dr) = SC/g, 

"Is = S/^m,P^ ^ ■ 

p Rewtonian viscosity, 

p density, [mL-3] 

T relaxation time, [t] 

<p friction factor 

<pQ = 64/Re 

cu angular velocity, 
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Terms in "brackets: 

U refers to plastic viscosity 

11 refers to apparent and local apparent viscosity 

Tig refers to S/Oj^^p 

^ refers to Nevtonian viscosity 

Subscripts : 

C capillary only 

e entrance 

m mean value 

P pipe or capillary 

p plug flow 

E rotational viscometer 

w wall of pipe or capillary 
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APPENDIX B 


NUMERICAL EXAMPLES 

The friction factors <p are to be determined for the flow of non- 
Newtonian materials in a smooth pipe line of_.__D = 2.0 centimeters and at 
different flow velocities. The flow velocities v are given in centi- 
meters per second. They are calculated from the specified required total 
mass -flow rates, the number of available pipe lines, the pipe diameter, 
and the density of. the materials, which is assumed to be 1.15 grams per 
cubic centimeter. 


Bingham Plastics 

The dashed curve in figure 9 is an assumed rotational- vis cometer 
flow curve of a Bingham plastic. The two flow parameters U and f 
are Independent of rate of shear. They are used to calculate Re[u] 
and PI, since Re[u] = Dvp/u and PI = fD/Uv. The friction factor <f> 
is then found from figure 3, The flow parameters are U = 0.19 poise 
and f = 400 dynes per square centimeter. Then the following values are 
obtained: 


V 

Re[u] 

PI 

<P 

Remarks 

q>o[u] 

50 

605 

84 

1.38 

Laminar 

0.105 

100 

1,210 

42 

.40 



.053 

200 

2,420 

21 

.12 



.030 

300 

3,630 

14.1 

.063 

1 


.040 

400 

4,850 

10.5 

.037 

Turbulent 

.039 

500 

6,050 

8.4 

.037 



.037 

1000 

12,000 

4.25 

.030 



.030 

2000 

24, 200 

2.1 

.024 



.024 


Pseudoplastic Materials 

The solid curve in figure 9 is an actual flow curve of a pseudo- 
plastic suspension measured on a concentric-cylinder rotational viscom- 
eter with R-b/Rc = 0.95. To obtain E, log is plotted against 

log in figure 10. Since no straight line results, N is not con- 

stant, but varies with rate of shear. The curve can be approximated by 
segments of linear functions, so that different constant values of N 
can be postulated for various mean-rate- of- shear intervals, as is indi- 
cated in figure 10. For more accurate values of N, the slope of the 
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logarithmic curve has to he determined at the point of specified local 
rate of shear. The local apparent viscosity is obtained from equation 
( 21 ). For extrapolations outside the flow-curve range of rate of shear,, 
the local apparent viscosity is obtained from figure 10 and equation (2l). 
It is used to calculate Re[ri], since Re[Ti] = Dvp/f]. The mean pipe rate 
of shear is Gj^^p/Cp,, where = 8 v/D. Then,, the local values of 

at which t] are to be obtained are equal to G^^pGp/Cp. The cor- 
relation factors C are obtained from figure 5 for the corresponding 
values of the local structure number N. The friction factor cp is 
found from figure 3 or equation (26). This procedure permits the calcu- 
lation of a.n local parameters. However, in case only the friction fac- 
tor is to be obtained, equation (27) can. be used, which eliminates all 
determinations of H and also the need for figure 10, unless extrapola- 
tion is necessary. The friction factor cp is again found from figure 3. 
The calculation of <p from equation (27) Instead of equation (26) has 
not only the advantage of greater simplicity but also of higher accuracy, 
since some uncertainties in N are eliminated. The calculated values 
are: 


V 


*local 

(a) 

Cp 

(t) 

%f/^p 

(f) 


Cor- 

rected 

®local 

(a) 

’'local 

(a) 

Rein] 

eq. ^e) 

(c) 

Remarks 

<PoLil ] 
(c) 

Rfihs] 

f, 

eq. (27) 

(a) 

50 

200 

4.00 

0.57 

350 

1.0 

350 

3.00 

o;9o 

128 

0.74 

Laminar 

0.49 

88 

0.72 

100 

400 

2.26 

.76 

525 

1.0 

525 

2.26 

.72 

320 

.26 



.20 

280 

.23 

200 

800 

2.26 

.76 

1060 

1.0 

1060 

1.90 

.49 

940 

.081 



.067 

815 

.078 

300 

1200 

2.26 

.76 

1580 

1.0 

1580 

1.66 

.40 

1,730 

.043 



.037 

1,480 

.044 

400 

1600 

1.66 

.86 

1860 

1.0 

I860 

1.66 

.37 

2,490 

.031 

TrauBition 

.030 

2,300 

.030 

500 

2000 

1.66 

.86 

2330 

1.0 

2330 

1.66 

,33 

3,500 

.040 

Turbulent 

.040 

3,220 

.040 

1000 

4000 

1,66 

.86 

4670 

1.0 

4670 

1.66 

.25 

9,250 

.033 



.033 

8,700 

.033 

2000 

8000 

1.66 

.86 

9350 

1.0 

9350 

1.66 

.18 

25,600 

.025 



,025 

23,000 

.025 


^■rom figure 10. 
^^From figure 5. 
<^001 figure 3. 
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TABLE I. - PIPE TRANSITION PRESSURE -LOSS COEFFICIENTS 
Cl for NEWTONIAN FLOW 


Abrupt contractions; 


D 2 /D 1 

0.8 

0.6 

0.4 

0.2 

0.0 

Cl* 

.13 

.28 

.38 

■ .45 

.50 (inlet from 
reservoir) 


Abrupt expansions; 


Cl = 



2 


Cl* = 1-0 ^outlet into reservoir, 




90° Bends: 


r/D 

1.0 

1.5 

2.0 

3.0 

4.0 

Cl* 

.40 

.32 

.27 

.22 

.20 


Valves (open) : 



Gate 

Globe 

Cl * 

0.2 

10.0 


90° Orifices (circular, slot ) : 


Cl = 



- 1 



0.0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 


Kfi 

.611 

.612 

.616 

.622 

.631 

.644 

.662 

.687 

.722 

.781 


fm 

00 

26,500 

1640 

319 

96.0 

37.5 

16.7 

7.9 

3.7 

1.5 

0 


Pipe entrance : 

Cl** = 2,16 (over entrance length Lg of pipe) 


*From ref. 1, pp. 156 and 49. 

**From ref. 10, p-. 48. 


















Rate of shear, dv/dr 
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Figure 1 . - Schematic flow curves of a Hewtonlan liquid. 
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Rate of shear; dv/dr 



Shearing stress; S '> — v — * Wall shearing stress; 

Plug flew 

Mean shearing stress; 

(a) Ideal. (h) Rotational viscometer. (c) Capillary viscometer 

and pipe line- 

Figure 2. - Schematic flov curves of a Bingham plastic. 
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?rl<jtlou factor, (from flov data and fig. 3) 
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Friction factor, ♦experimental direct preeiure meaaurementB) 


Figure 7 . - Friction f actora obtained for Blnglian. plaatlca from basic flov data and figure 3 com- 
pared vltb those from direct measure measurements. Pipe diameter, 0.75 to 3.05 inches; plastic 
Tiscoslty, 0.008 to 0.S35 poise; yield value, 0 to 276 dynes per square centimeter. 
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0 200 400 K)0 800 1000 1200 

Mean shearing stress, dynes/onip 


Figure 9. - Rotationai-viBcometer flow curves. 
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100 200 400 600 1000 2000 

Mean shearing stress, Sj^^jg, dynes/cm^ 


Figure 10. - Logarithmic plot of flow curve of 
pseudoplastic suspension (fig. 9) in accordance 
with eq-uation (19). 


NACA - Langley Field, Va 


